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1 Introduction 



Heavy ion collisions are an ideal tool to explore the QCD phase diagram. The goal is to study 
the equation of state (EOS) and to search for possible in-medium modifications of hadrons. By 
varying the collision energy a variety of regimes with their specific physics interest can be studied. 
At energies of a few GeV per nucleon, the regime where experiments were performed first at 
the Berkeley Bevalac and later at the Schwer-Ionen-Synchrotron (SIS) at GSI in Darmstadt, we 
study the equation of state of dense nuclear matter and try to identify in-medium modifications 
of hadrons. Towards higher energies, the regime of the Alternating Gradient Synchrotron (AGS) 
at the Brookhaven National Laboratory (BNL), the Super-Proton Synchrotron (SPS) at CERN, 
and the Relativistic Heavy Ion Collider (RHIC) at BNL, we expect to produce a new state of 
matter, the Quark-Gluon Plasma (QGP). The physics goal is to identify the QGP and to study its 
properties. 

By varying the energy, different forms of matter are produced. At low energies we study dense 
nuclear matter, similar to the type of matter neutron stars are made of. As the energy is increased 
the main constituents of the matter will change. Baryon excitations will become more prevalent 
(resonance matter). Eventually we produce deconfined partonic matter that is thought to be in 
the core of neutron stars and that existed in the early universe. 

At low energies a great variety of collective effects is observed and a rather good understanding 
of the particle production has been achieved, especially that of the most abundantly produced 
pions and kaons. Many observations can be interpreted as time-ordered emission of various particle 
species. It is possible to determine, albeit model dependent, the equation of state of nuclear matter. 
We also have seen indications, that the kaon mass, especially the mass of the K"*", might be modified 
by the medium created in heavy ion collisions. 

At AGS energies and above, emphasis shifts towards different aspects. Lattice QCD calcula- 
tions [U [2] predict the transition between a Quark-Gluon Plasma and a hadronic state at a critical 
temperature, Tc, of about 150 to 190 MeV at vanishing baryon density. The energy density at the 
transition point is about 1.0 GeV/fm^. It is generally assumed that chiral symmetry restoration 
happens simultaneously [1]. 

In the high-energy regime, especially at RHIC, a rich field of phenomena [3] has revealed it- 
self. Hot and dense matter with very strong collectivity has been created. There are indications 
that collectivity develops at the parton level, i.e. at a very early stage of the collision, when the 
constituents are partons rather than hadrons. Signs of pressure driven collective effects are our 
main tool for the study of the EOS. There are also strong indications that in the presence of a 
medium hadronization occurs through the process of quark coalescence and not through quark 
fragmentation, the process dominant for high-energy p-|-p reactions. 

We limit this report to the study of hadrons emitted in heavy ion reactions. The report is 
divided into two parts. The first part describes the phenomena observed from hadrons produced 
at low energies, whereas the second part concentrates on the search for signs of a partonic state at 
high energies. 

2 Hadron Production below 2 A GeV 

In the discussion of particle production in heavy ion collisions below 2 A GeV we focus on the 
most abundantly produced particles, pions and kaons. Pion production is the dominant channel 
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for particle emission, while kaon emission is a rare process. Kaons are produced during the high- 
density phase. Therefore, kaons are considered as ideal probes for the hot and dense fireball. 

In fixed-target p-|-p collisions the threshold for pion production is 0.29 GeV, for it is 1.58 
GeV and for K~ it is 2.5 GeV. The threshold for is lower than that for K~, because the energy 
balance has only to account for the excitation of a nucleon to a A in the case of K"*", whereas 
together with the K~ a has to be produced to conserve strangeness. 

The reaction dynamics of pions and kaons is quite different. The pion- nucleon cross section, 
shown in Figure [T] is large and thus pions are continuously absorbed through the A-resonance and 
re-emitted by its decay. This creation and disappearance can occur during the entire time evolution 
of the collision. The K'*"-nucleon interaction cross section, also shown in Figure [l] on the other 
hand, is small due to strangeness and energy conservation. There are no partners to react with 
and only elastic scattering and charge exchange can occur. Therefore, K"*" mesons are expected to 
leave the interaction zone early. These facts have led to the suggestion of time-ordered emission of 
these two species [Ij. The creation of K~ requires even higher energies. The dominant production 
mechanism is strangeness exchange. The interaction probability of K~ with nuclear matter is much 
higher than that of K"*" and the mean free path much lower. As a consequence, K~ cannot leave 
the reaction zone undisturbed and they are emitted later. 
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Figure 1: Elementary cross sections for K^p and vr^p interactions [5]. 

The properties of kaons are expected to be changed inside the nuclear medium [6l [7]. The 
scalar part of the KN potential of both kaon species is slightly attractive while the vector part acts 
in opposite direction. Thus, the K"*" N potential is slightly repulsive and the K~ N potential is 
strongly attractive. Potentials increase with density as shown in Figure |2] Kaons emitted in heavy 
ion collisions are an ideal tool to study this hypothesis. 

Kaons are also used to extract key parameters of the nuclear equation of state, like the nuclear 
compressibility. A large fraction of the work to extract the parameters of the EOS has been carried 
out through flow studies of the nucleons. A separate report in this Volume [8, describes these efforts 
in detail. 
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Figure 2: The effective mass of and of K~ resulting from the KN potential calculated within a 
relativistic mean field theory |7j as a function of the nuclear density [9] . 

Pion distributions were measured at the Bevalac in light symmetric systems {A < 40) and in 
asymmetric systems |10| [TT| [T2| [T3] and later also in the La+La system [14 . In these early mea- 
surements the slope of the high-py part of the spectra was taken as a measure for the freeze-out 
temperature, Tj^. This is a good approximation since in the blast-wave formalism [15] the radial 
expansion does not strongly modify the spectrum due to the low mass of the pions. The pion mul- 
tiplicity as a function of beam energy was measured by the Streamer Chamber Collaboration [16 . 
The difference between the measured pion yield and the yield predicted by cascade calculations pJT] 
has been taken as an early attempt to determine the equation of state of nuclear matter [18] . 

Positive kaons were measured early in the Bevalac program |19 | 120 1 1^. However, the statistics 
and systematics were not good enough to reconstruct a complete emission pattern the way it was 
done for pions. At Bevalac energies the production of negatively charged kaons is sub-threshold. 
The early measurements lack good statistics |i21j l22j , which adds to the problem to explain the 
observed phenomena close to the threshold where the production cross section varies rapidly with 
energy. Lambda and antiproton production [22| [23] IM] suffers from similar problems. A detailed 
overview of hadron production in high-energy nucleus-nucleus collisions at the Bevalac can be found 
in a review paper ^25j . 

Systematic, high-statistics measurements of pion and kaon production became feasible with the 
advent of the SIS accelerator at GSI. Two experiments, KaoS [26], a dedicated experiment for the 
measurement of K mesons, and FOPI [2^J, a large acceptance, multi-purpose detector, produced a 
solid body of beautiful results. References ^28^ and [29j give a detailed summary of pion and kaon 
production at SIS. 

Cascade codes are an ideal tool for the study of the mechanism of hadron production. At the 
Bevalac, pion production and kaon production [30^ was implemented in some models. Theory 
evolution not being part of this review, we limit ourselves to studying the production mechanism 
with the help (by means) of one particular model, the Isospin Quantum Molecular Dynamics 
(IQMD) model [HI [31]. Reference [32| reviews the theory development mainly concentrating on the 
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nuclear equation of state. We focus on reviewing multiplicities and spectral distributions. 



2.1 Multiplicities 

Pions are the most abundantly produced particles. They are created in individual nucleon-nucleon 
collisions via the A-resonance through the reaction A^A^ NA — > NNtt. Pions interact strongly 
with nuclear matter by forming baryonic resonances, e.g. vr+N— > A. The resonances themselves 
decay mainly by pion emission. Therefore, pions are expected to leave the collision zone over a 
long time span, mainly at a late stage of the collision when the system has expanded and cooled. 
For an experimental study of A-production we refer to Reference -VA . 

The pion multiplicity increases slowly as a function of the collision energy as shown in Figure [3] 
for inclusive collisions of a light (C+C) and a heavy (Au+Au) system. The ordinate M/A, the 
multiplicity per mass number A of one collision partner, accounts for the system size. The yield is 
higher for the lighter system. This reflects absorption of pions and holds only at these low energies. 
The kinematical conditions for tt"*" and tt~ are the same, the only difference is given by the N/Z 
ratio of the interacting nuclei. This effect is avoided here by plotting the sum of all three pion 
species. 




Figure 3: Multiplicity per mass number A (of one collision partner) for pions (diamonds, sum of 
all pion species), of (squares), and of K~ (circles) for Au+Au (filled symbols) and for C+C 
(open symbols) as a function of the beam energy. Values are taken from [33] . 

Figure [s] also shows the K+ and production yields for the two systems. As expected, the 
K"*" multiplicity is higher than the K~ multiplicity. The pion multiplicity is higher in the lighter 
system while for production the inverse observation is made. 

In central Au+Au collisions, densities of two to three times normal nuclear matter density may 
be reached [9l |3T1 132] . A sensitive probe to test such conditions is the production of strange hadrons 
at or below the production threshold. 
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Figure 4: Contribution of different production channels to the K"^ yield in central Au+Au collisions 
as a function of the beam energy [9] . 

The key mechanism for production close to the threshold is a multi-step process where the 
energy necessary for production is accumulated in intermediate resonances. Higher density increases 
the number of these collisions. Especially second generation collisions with sufficiently high relative 
momentum to create a occur most frequently during the high-density phase. Figure |4] shows 
the contribution of different production channels to the kaon yield in central Au-|-Au collisions 
according to IQMD calculations [9^. AN collisions dominate even above the corresponding NN 
threshold. 

The production process is quite different from K"*" and the production threshold is much 
higher. The strangeness-exchange reaction vrY K~N, with Y being A or S, represents an 
additional production possibility, as suggested by Ko [35] and demonstrated in References |36l |371 
[38] . It has a large cross section. The inverse channel causes the produced K~ to be absorbed |38] . 
As in the pion case, the succession of absorption and creation causes the K~ emission to be mainly 
in the late stage of the reaction. 

The 'K~ I'K^ ratio reflects the N jZ ratio of the colliding system. Assuming vr production via 
the A-resonance, the ratio is 1.95 (isobar value) for Au-t-Au collisions. The left panel of Figure [S] 
shows the vr^/vr''' ratio as a function of beam energy. The ratio increases with decreasing energy, 
exceeding even the isobar value. IQMD calculations [5] reproduce the trend but fail to describe 
the strong rise at very low incident energies. At very high energies the value reaches unity due to 
charge conservation [39] . 

The K^/K'*' ratio, shown in the right panel of Figure [s] rises strongly with increasing energy. 
This reflects the different production thresholds and kinematics for the two particles. The three 
systems studied here show quite similar values. The fact that different densities are reached in the 
different systems apparently does not matter. This trend is well described by transport models 
where strangeness exchange is the dominant channel and by statistical models [ID] . Again, at very 
high energies the value reaches unity [41] . 
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Figure 5: The tt~ /tt'^ ratio (left panel) and the K 
beam energy. The data are taken from References 
calculations 191. 



/K"*" ratio (right panel) as a function of the 
I28j and I29j. The sohd hne refers to IQMD 



In the following we study the multiplicities of pions and kaons as a function of the mass of the 
collision system, A + A, for inclusive reactions and as a function of the number of participants, 
Apart, for impact parameter selected reactions. 

Figure [g] shows the inclusive multiplicity per mass number, M/A, for pions as a function of A+yl 
at 1 A GeV incident energy. The pion multiplicity per A decreases strongly with A. This is in 
strong contrast to the behavior of production, shown in Figure [Tj At 1 ^ GeV, the normalized 
kaon yield rises by a factor of about 3. Figure [t] also shows the variation of multiplicities of K"*" 
mesons per A with incident energy as well as those of K~ mesons at 1.5 A GeV. The lines are 
functions M ~ A"' fitted to the data with the resulting values for 7 given in the figure. For K"*" 
production the extracted values of 7 decrease with incident energy. This reflects the decreasing 
influence of intermediate energy storage via A as less multiple collisions are needed. Considering 
the much higher threshold for K~ production, one would expect at the same incident energy a 
much stronger rise. However, at 1.5 A GeV the values of 7 for and for K~ production are 
about equal, demonstrating that production and production are strongly correlated. 

Figure [s] shows M/ Apart for Ni+Ni and Au+Au collisions at 1.5 A GeV as a function of Apart- 
It demonstrates that the multiplicities of both kaon species exhibit the same rise with the number 
of participating nucleons despite the fact that the thresholds for the production of the two particle 
species are very different. This observation again confirms that K"*" and productions are corre- 
lated. As is also shown in Figure [8] the pion multiplicity per Apart as a function of Apart is rather 
fiat. This observation has been made at several incident energies |16| H5] . 

It is important to distinguish the trends observed in the quantity M/A for inclusive reactions 
and M/Apart studied for different centralities. While for kaons these trends are quite similar as 
demonstrated in Figures [7] and [Sj it is not the case for pions. While M/Apart is flat, M/A decreases 
when comparing different systems (see Figure [g]). Consequently, M/Apart for different systems 
gives different values. Results from inclusive and from centrality-selected studies can not easily be 
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Figure 6: Pion (sum of vr^, vr'', vr") multiplicities per A for inclusive collisions of p+p (divided by 
2), C+C, Ni+Ni and Au+Au reactions as a function A + A, the total mass of the system. The 
data are taken from References |28t W2\ US] . 
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Figure 7: Multiplicities per mass number Mj A as a function of A for K"*" (full symbols) and for K 
(open symbols) for inclusive C+C, Ni+Ni, and Au+Au reactions. The lines represent the function 
M ^ A^ fitted to the data. From Ei . 
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Figure 8: Dependence of the multiplicities of K^" (upper panel) and of K~ mesons (middle panel) 
as well as of pions (lower panel) on Apart- Full symbols denote Au+Au, open symbols Ni+Ni, both 
at 1.5 ^4 GeV . The lines are functions M ~ A^^^^ fitted to the data separately for Au+Au (solid 
lines) and Ni+Ni (dashed lines). The data have been measured at ^lab = 40°. Taken from [29j. 
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compared. 

The similar rise of both, and K~ , as a function of the collision centrality shown in Figure [s] 
and as a function of the system size in Figure [7] suggests that the production mechanisms of the 
two kaon species is correlated. The K~ in heavy ion collisions at SIS energies are mainly produced 
via strangeness-exchange. On the other hand, strangeness has to be conserved when producing 
these hyperons and the energetically most favorable way is to produce them together with K"*" (or 
mesons. Thus the production of and of K mesons is coupled via the strangeness-exchange 
reaction and the K~ inherits the same dependence on the system size and on the collision centrality. 
In Reference ^K] it is argued that the strangeness-exchange channel reaches chemical equilibrium. 
Thus the K^/K+ ratio is proportional to the pion density, both at SIS and AGS energies. 

As a consequence of Figure [s] the ratio K~/K+ as a function of Apart is constant and it is 
the same for Au-|-Au and for Ni-|-Ni. Transport models having strangeness exchange as dominant 
production channel for K~, describe this trend fairly well [38]. The statistical model reproduces 
the measured values using nominal masses for all particles [ID]. This model predicts a value 
independent of centrality and of system size. At these low incident energies, a canonical treatment is 
required leading to multiplicities of strange particles depending on the size of the fireball (correlation 
volume). In this approach the size dependence of both, and K~, is the same and a constant 
K~/K+ ratio results. 



2.2 Spectra 

The spectral distributions of the different particles in general contain information about the freeze- 
out condition and about the collective expansion of the system created in heavy ion collisions. 

The invariant cross sections for pion production in three different systems are shown in Figure [9] 
as a function of the center-of-mass kinetic energy. Above 0.4 GeV kinetic energy in the cm. 
system all spectra can be described by a Boltzmann distribution (exponential function in the given 
presentation). The spectra also exhibit an enhancement at the low-energy end if an exponential 
function is fitted in the high-enery part. This excess is stronger for heavy systems. There are 
different possible explanations for this enhancement. Late time pion emission from A-resonance 
decay |13j and continuous pion emission [l8] have been proposed. The most detailed investigation 
to date has been done by Weinhold et al. [l9] . 

As mentioned earlier, the ratio of the tt~^ and vr^ multiplicities reflect the N/Z ratio of the 



colliding nuclei and the opposite Coulomb force. Figure 10 shows the double differential cross 
sections of negatively and positively charged pions from central collisions of Au-|-Au nuclei at an 
incident beam energy of 1 ^4 GeV. The yield of vr" is higher due to the neutron excess. The different 
shapes demonstrate the boost caused by the Coulomb force. Similar observations have been made 



at the Bevalac when studying forward emission of pions |50j . The right panel of Figure 10 compares 
the 'K~ jix^ ratio as a function of the kinetic energy of the pions with calculations using a static 
Coulomb potential [U] which is not able to describe the dynamics of pion emission. For pions with 
lower kinetic energy a weaker Coulomb force seems to act. This can be understood by the time 
sequence of pion emission; high-energy pions are emitted earlier than low-energy pions. 



Both and K~ are produced in a complicated sequence of interactions. Figure 11 shows 
K"^ and K~ spectra at mid- rapidity as a function of the kinetic energy i?c.m. — m^c? for three 
different systems and various beam energies [29j. The spectra have a Boltzmann shape to a very 
good approximation. The inverse slope parameters of the are always higher than those of the 
K~ |5H l52l [53] [29] . Heavy systems exhibit shallower inverse slopes than lighter ones. 



10 



■« 10' 

"> 

O 

"a 

s 1 

a 

'"■o 10 " 



71"*" Au+Au 



It* Ni+Ni 



;p'. 'l.SAGeV 

□ 1.0 AGeV □ 1.0 AGeV 

■O.SAGeV 



central 



central 



It* C+C 




10' 

o 

a 

'ts 10 " 









° 7i" Au+Au 


It" Ni+Ni 


It" C+C 




v. • 1.8 AGeV 


A • 2.0 AGeV 


% ° 1.0 AGeV 

\ 
% 

\ 


\=X = 1.0 AGeV 


% = 1.0 AGeV 

V 


\ 






central 


central 


inclusive 









C„(GeV) 



Figure 9: Spectra of positively (upper panel) and negatively (lower panel) charged pions in the 
center-of-mass frame in a Boltzmann representation for various reactions |47| H3] . 
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Figure 10: Left panel: Double differential cross sections of negatively and positively charged pions 
from central collisions of the reaction Au+Au at an incident beam energy of 1 ^4 GeV and at an 
emission angle of ^lab = 44° ±4° . Right panel: Ratio of vr^/vr^ as a function of the kinetic energy 
of pions for Au+Au collisions. From j47j 
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Figure 11: Inclusive invariant cross sections as a function of the kinetic energy Ec.m. —mQC^ for 
(left panel) and for K~ (right panel) for three systems and various beam energies at mid-rapidity 
(6»c.m. = 90° ±10°). From [29]. 
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Figure 12: Influence of rescattering of the K"*" and of the repulsive KN potential for central Au+Au 
collisions at 1.5 A GeV. Left panel: Initial and final distributions of kaons which have scattered 
twice or more. Right panel: Influence of the KN potential on the spectral shape demonstrated by 
selecting kaons that never scattered {Nq = 0) and comparing the initial and final spectra [9]. 



It is interesting to study the different effects contributing to the K"*" slope in an IQMD study. 
Initially, the energy available for production is low and the slope of the spectrum at creation 



is therefore quite steep. Figure 12 shows the effect of re-scattering and of a possible kaon potential 
on the final slope of K^. Re-scattering makes the slope shallower. This effect increases as the 
beam energy increases and as the mass of the system increases as there is more scattering in 
heavier systems. The predicted repulsive K"*" potential would increase the production threshold. 
When leaving the system, the K"*" has to acquire its nominal mass and the excess energy serves 
to accelerate it. This effect, however, is small and only visible at low momenta. For experimental 
reasons low-momentum kaons can be studied best with neutral K'^ mesons. The corresponding 
experiments are still ongoing |54| |55] . 

Re-scattering has a strong influence on the K"*" spectra. The spectral shape of the K~ mesons 
is also influenced by scattering. However, K~ are mainly absorbed when interacting with nuclei 
and absorption will change the slope in a momentum-dependent way. The low-momentum K~ 
undergo stronger absorption which leads to an increase in the apparent temperature. The various 
influences on the shape, resulting from re-scattering, absorption, and the attractive KN potential 



are demonstrated in Figure 13 

Figure 14 gives a compilation of the measured inverse slopes of inclusive vr^ (only high-energy 
part) and inclusive K"*" for C-|-C and Au-|-Au collisions as a function of collision energy, clearly 
showing the trends discussed before. The inverse slopes increase with beam energy and at the 
same energy those of the heavier systems are higher. It is interesting to note that the values of the 
slopes of the high-energy pions and of the K"*" agree. They seem to reflect the temperature of the 
environment, if this term can be used at all as a global equilibrium is not achieved. Those of K~ 



deviate as already shown in Figure 11 



Finally, we investigate the impact parameter dependence of the kaon spectra. Figure 15 (left 



panel) shows the spectra measured by the KaoS Collaboration for Au-|-Au reactions at 1.5 ^4 GeV 
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Figure 13: Left panel: Influence of absorption only demonstrated for Au+Au at 1.5 A GeV 
based on IQMD calculations. Right panel: Effect of scattering of K~ and of the attractive KN 
potential |9]. 
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Figure 14: Measured inverse slopes of high-energy pions (left panel) and of mesons (right panel) 
from inclusive C+C and Au+Au collisions at mid-rapidity |29 | 134 ) . 
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Figure 15: Left panel: Invariant cross sections for Au+Au at 1.5 A GeV close to mid- rapidity 
for different centralities. The uppermost spectra correspond to the most central collisions. The 
subsequent bins are shown from top to bottom with decreasing centrality. The lines represent fits 
to the KaoS data [29 . Right panel: Summary of inverse slope parameters [29] . 
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Figure 16: Left panel: Time profiles for production (dashed lines) and last contact (solid lines) 
of K"*' (top) and K~ (bottom) mesons. Only those K~ are shown which finally leave the system. 
Right panel: Density distribution at the point of production (dashed lines) and at the point of 
last contact (solid lines) for K"*" (top) and K~ (bottom) mesons. The simulations are for central 
Au+Au collisions at 1.5 A GeV. Figures are taken from [9]. 



for different centralities for K+ and K^. The inverse slope parameters, T, are evaluated as a 
function of centrality or Apart- This demonstrates that the slope parameters are weakly depending 
on centrality and that the slope of K~ is always about 20 MeV lower than that of K"*". 



2.3 Conclusions 

Heavy ion collisions at 1 to 2 A GeV incident energy probe nuclear matter at about two to three 
times its normal density. The production and emission of K"*" andofK" exhibits distinct differences. 
The key points are summarized in a model study showing the emission time of the two kaon species 
and the corresponding density profiles of the medium within the IQMD approach [9 . The left panel 



of Figure 16 shows the time distributions at production and at last interaction. K are produced 



later than K"^ and they also leave the system later. The right panel of Figure 16 shows the density 
of the medium at production and at last contact. The bulk of the K"*' is produced when the 
density is twice nuclear matter density. K"*" are emitted predominantly at a density of p = 1.5 po- 
Their yield cannot be changed between production and emission due to strangeness conservation. 
Therefore, the K"*" yield carries information from the high-density phase. As K"*" production is at 
or below the NN threshold, its yield is very sensitive to the dynamics in the fireball. K"*" production 
serves as an ideal tool to extract information on the nuclear EOS. In order to avoid uncertainties 
both experimentally and theoretically, it is useful to study ratios of multiplicities [Mj, e.g. the 
ratio of multiplicities from Au+Au and C+C collisions. From this ratio a rather soft EOS with a 
compression modulus of around 200 MeV has been extracted [56 1 [57] . This is in agreement with 
results from flow studies [HI 1^ . 

The observed K~ mesons are predominantly produced and emitted late and from a region 
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Figure 17: Mid-rapidity (\y\ < 0.5) proton (left column) and antiproton (right column) transverse 
mass distributions from Au+Au collisions at y^s^ = 130 GeV. From bottom to top the curves are 
ordered according to centrahty bins, 70-80%, 60-70%, 50-60%, 40-50%, 30-40%, 20-30%, 10-20%, 
and 0-10%. Results from thermal model fits are shown as dashed lines. 



of density below nuclear matter density, whereas the K"^ are produced at twice nuclear matter 
density. It is obvious that heavy ion reactions in this energy range are well suited to study 
potential effects. This is not the case for K~ production, since the effects are expected to be 
small. Other observables, like the azimuthal emission pattern, also show sensitivity to in-medium 
interactions [591 EO] . 



3 Hadron Production from AGS to RHIC 

In this section we discuss the evolution of spectra of identified hadrons as the center-of-mass energy 
is increased from AGS to RHIC. For AGS and SPS we will summarize the spectra in terms of 
inverse slope parameters, T, and show how radial flow develops as a function of energy. At RHIC 
energies we will use the wealth of precise data on radial and elliptic flow to argue that collectivity 
develops at the partonic level. 

We start with a discussion of the transverse momentum spectra and the way they are ana- 
lyzed and important dynamical parameters are extracted using the example of spectra from RHIC 
experiments. The methods developed are applicable to all energies discussed here. 



Figure 17 shows mid-rapidity proton and antiproton transverse mass distributions for several 
centrality bins. The data are from Au+Au collisions at ^/s^ = 130 GeV [61^. The transverse 

mass, rriT, is given by mx = \/{p'^ + mQ), with tuq the rest mass of the proton (antiproton). 



It is evident that the distributions become more convex as centrality increases. The change in 
shape is an indication for strong radial flow. In the presence of collective flow the transverse 
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Figure 18: Mid-rapidity transverse mass distributions from Au+Au and p+p collisions at ^/s^ 
= 200 GeV for mesons (left panel). Dashed lines represent exponential fits to the data and the 
dot-dashed line is the power-law fit to the p-|-p spectrum. In the right panel H particles are shown. 
Dashed lines represent exponential fits. The data are from STAR [66] . 



mass distributions, especially for heavy mass particles, will not have a simple exponential shape 
at low transverse mass [15]. This effect becomes particularly strong when the temperature is low 
compared to the collective velocity. Two-parameter fits |62j of the freeze-out temperature, Tf^, 
and the collective velocity, f5t, motivated by hydrodynamics, are useful to separate the collective 
motion from random thermal motion in the measured spectra. Here we apply blast wave fits 
from Reference [63] . These fits simultaneously describe experimental spectra of charged pions [M] , 
kaons [65^, protons and antiprotons. A velocity profile (3t{i") = Ps{f/R)^'^ is used in the fit, where 
R and (3s are the radius and the surface velocity of the source, respectively. The fit-results are 



shown as dashed lines in Figure 17 The collective velocity parameter < f3t > increases from about 



0.41 c to 0.55 c for the most central collisions. 



Figure 18 shows </) meson and H hyperon transverse mass distributions from 200 GeV Au-|-Au 
and p-|-p collisions. The spectral shapes are quite different from the ones shown in Figure [TT] All 
spectra with the exception of the (j) niesons emitted in p-|-p collisions have an exponential form: 

fexp = A ■ e , (1) 

where A and T are the normalization constant and inverse slope parameter, respectively. The 
dashed lines in the figure represent exponential fits. In the case of (p rnesons emitted in p-|-p 
collisions, a power-law fit gives better results: 

/p = ^-(1+Pt/po)-" (2) 
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Figure 19: Inverse slope parameter, T, of pions, kaons, protons, and deuterons from central collisions 
at the AGS {Js^ ~ 5 GeV) for p+Au , Si+Al , Au+Au . 



Here A is the normalization constant and po and n are free parameters that describe the shape 
of the distribution. The dash-dotted line represents the fit. In general, mesons emitted from p+p 
collisions are following a power-law type distribution with a pQ of about 12 GeV/c and a n of six, 
while the baryons are closer to an exponential function [ 67] . 

3.1 Systematics of Spectra and Slopes at AGS and SPS 

The AGS program has yielded a large body of emission spectra mainly for pions, kaons, and protons. 
At the SPS there is an even larger wealth of data. Rather than reviewing the data in detail, we will 
concentrate on summarizing the emission spectra in terms of inverse slope parameters or apparent 
temperatures. 



Figure 19 shows the systematics of the extracted inverse slope parameters as a function of the 



mass of the emitted particles at AGS energies for three systems with different total mass. The 



data are taken from the E802 experiment |6§J . Figure 20 gives the corresponding systematics 
for SPS energies. Here the data are taken from the NA44 experiment [62]. At both energies a 
clear trend emerges: a linear relationship between inverse slope parameter and particle mass. This 
dependence can be parameterized as a freeze-out temperature plus a linear term in particle mass 
with a coefficient that is related to the flow velocity. The p+p and the p+nucleus systems do not 
show any dependence on mass. Thus, in this approach, they do not show collective expansion. For 
the heavier systems we observe a collective expansion velocity that increases with the mass of the 
system. A freeze-out temperature of 140 to 150 MeV has been extracted. The systematics observed 



in Figure 19 and in Figure EO^ has been taken as a strong indication of the existence of collective 



radial expansion in heavy ion collisions in this energy range |691 [70] . 



Figure 21 shows in more detail again the inverse slope parameter as a function of particle mass 



for central Pb+Pb collisions at 158 GeV with data from NA44 [62], NA49 [71], and WA97 [72]. The 
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Figure 20: Inverse slope parameter, T, of pions, kaons, protons, A, and H from central S+S and 
Pb+Pb collisions and minimum bias p+p collisions at the SPS {^/s^ ~ 17 — 20 GeV). 
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Figure 21: Measured inverse slope parameter, T, (filled triangles) and RQMD model predictions 
(open circles) as a function of particle mass of central Pb+Pb collisions at ^/s^ ~ 17.3 GeV. The 
line serves to guide the eye. 
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linear increase from pions to deuterons is indicated by the solid line. The inverse slope parameters 
for multi-strange baryons and for (p mesons clearly deviate from this general trend. Those particles 
do not acquire the flow velocity appropriate to their mass during the expansion. They interact 
weakly with the expanding system since their hadronic cross sections are reduced compared to 
non-strange hadrons . The observed trend is confirmed by RQMD calculations [71] , also shown 



in Figure 21 In RQMD, the interactions are modelled by hadronic scattering alone. The consistency 
between the model and experimental data clearly indicates that this trend is due to scattering in 
the hadronic phase. 

For light hadrons, the bigger their mass the higher the value of the inverse slope parameter. 
The multi-strange hadrons like Xi and Omega, on the other hand, do not follow this trend. Their 
inverse slope appears to decrease. As discussed in Reference [73], this different behavior may be 
caused by the small hadronic cross sections of the multi-strange hadrons. They decouple early from 
the system and do not acquire the degree of collectivity characteristic for the abundantly produced 
light hadrons. As we will see later, the situation will change when the beam energy is increased at 
RHIC. 

3.2 Spectra and Radial Flow at RHIC 

In the previous Section, we compared the systematic results from AGS and SPS with predictions 
from transport models. In case of the transport approach, the driving force is rescattering. No 
assumption is made regarding the equation of state and the nature of thermalization. A sufficient 
number of rescatterings will naturally lead to collectivity and eventually to local equilibrium. We 
can also look at the collisions from the point of view of hydrodynamics. In Figure [22] we compare 
the data with hydrodynamic model calculations [75]. In the following, we will assume that we can 
identify bulk properties by measuring the most abundantly produced particles, like pions, kaons, 
and protons. 

The values of the mean transverse momentum of bulk particles, from PHENIX [76] (triangles) 



and STAR [3^ are shown in the left panel of Figure 22 as a function of collision centrality. Also shown 



in the plot are the hydrodynamic model calculations where the bands indicate the uncertainties 
caused by the initial conditions. Within error bars, the data for bulk particles are well reproduced 
by the model calculations, implying the validity of the application of hydrodynamic approach for 



bulk productions in high-energy nuclear collisions. In the right panel of Figure 22 we compare the 
< > of hadrons as a function of hadron mass for pions, kaons, protons. A, H, and 0, from central 
Au-|-Au collisions. In the mass region below 1 GeV the hydrodynamic calculations can account 
for the measured results reasonably well. In the mass region above 1 GeV the measured mean 
transverse momenta are below model predictions. In this region most of the produced hadrons 
contain one or more strange quarks. The interaction cross sections of strange (or multi-strange) 
hadrons are smaller than those of non-strange hadrons j73j . As a result, they decouple from the 
system relatively early. Therefore the values of the < pT > are lower than those from hydrodynamic 
model calculations which can reproduce the mean transverse momentum for bulk particles, but 
over-predict < px > for multi-strange hadrons. 



Figure 23 shows contours in the temperature versus collective velocity plane. Dashed and 
solid lines represent 1-a and 2-a contours, respectively, extracted from blast wave fits [63] . Numbers 
on the top indicate the centrality selection. For bulk particles 9 centrality bins are presented. For (p 
and 0, only the results from the most central collisions are shown. Contours from p-|-p collisions [77j 
are also shown. 
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Figure 22: Average transverse momenta for pions, kaons, and protons as a function of collision 
centrality (left panel) and from central collisions as a function of the hadron mass (right panel). 
Hydrodynamic model results ^75j as shown as bands. 
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Figure 23: contours from blast wave fits for vr, K, and protons and for the multi-strange hadrons, 
<j) and O. Numbers on top indicate the centrality selection. Dashed and solid lines represent 1-a 
and 2-0" contours, respectively. 
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Figure 24: Inverse slope parameter, T, as a function of hadron mass for Pb+Pb central collisions 
at SPS (open circles) and Au+Au central collisions at RHIC (filled circles). The hatched area 
indicates the expectations discussed in the text. 



With increasing centrality the bulk temperature decreases and the collective velocity increases. 
The velocity becomes as high as 60% of the speed of light for the most central collisions. In contrast, 
the multi-strange hadrons, cp and O, emitted in central collisions have a temperature of T ~ 180 
MeV and an average velocity of /3 ~ 0.4. From Figure 18 we also know that their temperature and 
expansion velocity is not sensitive to the collision centrality. 

The thermal freeze-out temperature for the bulk is about 100 MeV. For the multi-strange 
hadrons it is about 170 MeV, the same as the bulk chemical freeze-out temperature [701 ES ESI 
[HO] IHT] . This temperature is close to the value of the expected phase transition temperature [2] . 
Again, we can conclude that the multi-strange hadrons do not participate in the full development 
of the expansion, they decouple from the system near the hadronization point, T ~ 170 MeV and 
/3 ~ 0.4. A mechanism that could lead to early freeze-out is the effect that the hadronic cross 
sections for strange and especially multi-strange hadrons are reduced [73], as explained earlier. 

An alternative scenario has been proposed by Broniowski and Florkowski [82]. They assume 
a single freeze-out and explain the apparently low freeze-out temperature and the large collective 
velocity for the bulk as an effect of resonance decay. However, it is not clear that this model can 
also explain the data on resonance production |j83j . 



Figure 24 shows the inverse slope parameters of identified hadrons as a function of particle mass. 
At SPS energies the multi-strange and charm hadrons do not show an increase as mass increases. 
This is indicated by the solid band. At RHIC energy, however, the inverse slope parameters for 
H and Q seem to increase with mass, indicating that the collectivity developed already in the 
partonic phase [7D]. The hatched band indicates where future data points would have to fall for 
this assumption to be correct. 
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Figure 25: Left panel: The V2{pt) for tt, Kq, p, and A from minimum bias Au+Au collisions at 
y/s^ = 200 GeV [85]. Right panel: V2 and px scaled by the number of constituent quarks. The 
dashed lines represent the elliptic flow of vr, K, p, A, H, and J7 from hydrodynamic calculations. 
Dot-dashed lines are the fits to Kq and A V2 distributions. 



3.3 Partonic Collectivity at RHIC 

More conclusive evidence for flow at the parton level can be derived from elliptic flow measurements. 
Elliptic flow is described in a separate article of this Volume pi] . We will only present a limited 
set of data here. 



The left panel of Figure 25 shows the measured elliptic flow, V2, from Au+Au collisions at ^/s^ 
= 200 GeV for tt, K^, p, and A from PHENIX [85- and STAR [Sj. From top to bottom the dashed 
lines represent the elliptic flow of vr, K, p, A, H, and from hydrodynamic calculations [50] • At 
low pt, the hydrodynamic calculations can well reproduce the V2 measurements. At higher pT, V2 
becomes saturated and hydrodynamic results over-predict the data. The baryons saturate above 3 
GeV/c with V2 ~ 0.2. The mesons saturate at lower values of V2- When scaling the values of V2 
and pt with the number of constituent quarks (NCQ) of the corresponding hadrons, all particles 



should fall on a single curve, as demonstrated in the right panel of Figure 25 , This was predicted 
by coalescence models [57] . 

Scaling with the number of constituent quarks or with the number of valence quarks has been 
taken as strong empirical evidence that the constituents that acquire collectivity in the expansion 
are partons (quarks) and not hadrons. 



The V2 values in Figure 25 are fitted with an empirical functional form. The fit results for 
and A are shown as dot-dashed lines in both panels of Figure [25} For kaon, proton, and lambda 
the scaling works well for small transverse momenta {px/nq < 2.5 GeV/c). Pions do not follow 



25 





^NN (GeV) 



Vs" M (GeV) 



Figure 26: Energy dependence of the elliptic flow V2 from minimum biased Au+Au or Pb+Pb 
interactions (left panel) and of the radial flow velocity parameter {Pt) from central Au+Au or 
Pb+Pb collisions (right panel). 



the scaling. A large fraction of pions from heavy systems is not produced directly, but comes from 
the decay of resonances |83t I88j . At mid-rapidity as much as 80% of pions are from resonance 
decays. The dominant sources for pion production are p, K*, and baryon resonances like A. 
When taking resonance decays into account, the pion V2 follows the scaling with the number of 
constituent quarks |89j. 

The scaling of the transverse momentum dependence of V2 with the number of constituent 
quarks of the produced hadrons is a powerful indication that the constituents that acquire flow are 
not the hadrons, but partons at an early stage of the collision. This is the strongest indication so 
far for deconfinement and a hint for the existence of a Quark-Gluon Plasma. 

Let us close this section by showing the energy dependence of the collective flow observables 



in Figure 26 The elliptic flow V2, and the radial flow parameters, (5t, increase as a function of 
collision energy. Above 6 GeV f 2 increases monotonically implying that more and more early flow 
develops as the collision energy increases. Assuming that hadronic rescattering effects have already 
been maximized at SPS energies (y^s^ ~ 17 GeV), the net increase in the transverse velocity at 
RHIC should be due to partonic interactions. 

3.4 Outlook 

In the near future, it is important to quantify the partonic collective flow with high statistics data of 
V2 measurements for all hadrons. Elliptic flow of (j) mesons will be important because (p mesons are 
produced through quark coalescence and not through the kaon fusion channel [90]. These results 
will provide direct information on the partonic phase. 
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Gauging thermalization is an important next step in the quest for the equation of state and 
for the existence of a Quark-Gluon Plasma. In lattice QCD the QGP temperature at top RHIC 
energy is of the order of 0.3 to 0.5 GeV |9H l92] . The mass of a charm quark is much heavier than 
this temperature. This means thermal production of charm quarks is negligible. Charm quarks are 
produced in first chance collisions. Thermalization of the light quarks might be tested by measuring 
the elliptic flow of heavy quarks. Heavy quarks can develop flow only if they are dragged along 
during the expansion by frequent collisions with light quarks Frequent interactions amongst 

light quarks will lead to thermalization. The measurement of a finite V2 value for electrons from 
the semi-leptonic decay of heavy-flavor mesons at high transverse momentum [9l] cannot yet be 
taken as indication of charm quark flow in the hydrodynamic sense. 

It is also very important to pursue a temperature measurement with thermal photons and di- 
leptons. NA60 |95j recently obtained "temperatures" for di-muons from the continuum that fit 



well in the mass systematics of Figure 24 , Thermalization of matter is a necessary condition for 
measuring the temperature of the QGP. The temperature is yet another unknown parameter in the 
determination of the equation of state. Thermalization is an important condition for our ability 
to explore the QCD phase diagram, for example by searching for the tricritical point [96] in the 
upcoming energy scan program at RHIC and later at FAIR. 



4 Summary 

In heavy ion collisions at energies of 1 to 2 A GeV densities of two to three times normal nuclear 
matter density are reached. K"*" mesons production is of special interest as K+ are produced close 
to the NN threshold. Their production yield is very sensitive to the energy balance in the fireball 
and thus to the properties of the nuclear equation of state. K"*" interact only elastically and their 
yield carries information from the early, high-density stage of the collision. From the comparison 
of the measured yields with model calculations the stiffness of nuclear matter has been extracted. 
Pions and interact strongly with the medium. They are emitted late in the expansion phase. 

Kaon emission close to the threshold is a sensitive tool to identify effects of in-medium KN 
potentials. Potentials might cause a change of the effective thresholds and of the yields. Presently, 
experimental data and the understanding of the model parameters are not good enough to distin- 
guish the effects of potentials from other effects, like scattering and absorption. and K*^ spectra 
at very low momenta are expected to be more sensitive. With the upcoming FAIR facility at GSI 
it will be possible to study potential effects through high statistics measurements by including also 
charm production. 

At the higher energies strong, pressure driven collective effects have been observed. The fact 
that elliptic flow scales with the number of constituent quarks and that the multi-strange hadrons 
decouple from the system early have been taken as strong indication that collectivity develops in 
the partonic stage and that hadronic interactions play a minimal role for the development of elliptic 
flow. We can conclude that a partonic state has been created. 

Deconfinement is not sufficient for the existence of a QGP. Equally important is the concept 
of thermalization. The heavy flavor program at RHIC and at the LHC will shed new light on this 
important question. 
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